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Using the model independent method of Villante, Fiorentini, Lisi, Fogli, Palazzo, and the rates 
measured in the SNO and Super-Kamiokande solar neutrino experiment, we calculate the amount 
of active v^ or v T present in the flux of solar neutrinos on Earth. We show that the probability 
of i/ e — » Vy,, T transitions is larger than zero at 99.89% CL. We find that the averaged flux of ^ M , r 
on Earth is larger than 0.17 times the 8 B v e flux predicted by the BP2000 Standard Solar Model 
at 99% CL. We discuss also the consequences of possible v e — > z/ M]T or v e — > v a transitions of solar 
neutrinos. We derive a model-independent lower limit of 0.52 at 99% CL for the ratio of the 8 B z/ e 
flux produced in the Sun and its value in the BP2000 Standard Solar Model. 
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The first results of the SNO solar neutrino experiment Q have beautifully confirmed the existence of the solar 
neutrino problem. A comparison of the neutrino flux measured through charged-current interactions in the SNO 
experiment with the flux measured through elastic scattering interactions in the Super-Kamiokande experiment 0] 
f»0 ■ shows an evidence of the presence of active v^ or 1 v T in the solar neutrino flux measured by the Super-Kamiokande 
^ | experiment p], pj. Such a presence represents a very interesting indication in favor of neutrino physics beyond the 
Standard Model, most likely neutrino mixing that generates oscillations between different flavors (see B). 

The purpose of this paper is to quantify the amount of this flux of active v^ or v T in a model-independent way in 
the framework of Frequentist Statistics 2 . 

The authors of Refs. flq, pj have noted that the response functions of the SNO and Super-Kamiokande (SK) ex- 
periments to solar neutrinos can be made approximately equal with a proper choice of the energy thresholds of the 
detected electrons. It turns out that given the threshold T^ NO — 6.75 MeV, the two response functions are approx- 
imately equal for TJ? K — 8.60 MeV ||. In this case the SNO and Super-Kamiokande event rates normalized to the 
BP2000 Standard Solar Model (SSM) prediction M can be written in a model-independent way as B 



Ph! RsNO = fn(P^» e ), (1) 

*f-j ' ^sk = /b <P„.^,.) + /b V^V < P ".-«V.r> . (2) 

> ; ^'t 

where /b is the ratio of the 8 B v e flux produced in the Sun and its value in the SSM ||, (P Ve -> V e) is the survival 
Vh ' probability of solar v e J s averaged over the common SNO and Super-Kamiokande response functions, 



( CT ^,.) 



0.152 (3) 



is the ratio of the averaged u^. T and v e cross sections in the Super-Kamiokande experiment, and (P Ve ^ v T ) is the 
averaged probability of v e — > i/^ iT transitions. 
Calling 

Ra = Rsk — Rsno > (4) 



'Electronic address: 5iunti@t0.infn.it; URL: tittp://www. to. infn.it/~giunti 

1 In this paper the conjunction "or" is used as a logical inclusive disjunction (the sentence is true when either or both of its constituent 
propositions are true). 

2 Since the results that we obtain are not too close to physical boundaries for the quantities under discussion and we assume a normal 
distribution for the errors, the numerical values in the framework o£_Bayesian Probability Theory with a flat prior are close to those 
obtained here, but their meaning is different (see, for example, Ref. [B|). 



from Eqs.(||) and (||) we have 



RA = fB^4{Pu.^, T ). (5) 



Therefore, Ra is the rate of ^ MjT -induced events in the Super-Kamiokande experiment, relative to the i/ e -induced rate 
predicted by the SSM. 

Considering the data of the Super-Kamiokande experiment above the energy threshold T e SK = 8.60 MeV and the 
BP2000 Standard Solar Model H, the measured values of i?sNO an d -Rsk are: 

R%* = 0.347 ±0.029 [1], (6) 

R c s x i = 0.451 ± 0.017 §,§. (7) 

Adding in quadrature the uncertainties of i?sNO and .Rsk, for Ra we obtain: 

? oxp 



i^ xp = 0.104 ±0.034. (8) 



The standard deviation of R°? p is 



and we have 



cxp 



0.034, (9) 



Ej>CXp 

-^r = 3.06±l. (10) 

a A 

Hence, the central value of Ra is 3.06cr away from zero, implying an evidence of solar v e — > ^ M/r transitions [EL p). 
Our purpose is to quantify the probability of these transitions and possibly derive a lower limit. 

The authors of Ref. \M calculate the probability of a fluctuation larger than the observed one assuming Ra = 0: 
for normally distributed errors the probability of a fluctuation larger than 3.06er from the mean is 0.11%. 

Recently some frequentist methods have been proposed that allow to obtain always meaningful confidence intervals 
with correct coverage for quantities like Ra that are bound to be positive by definition ^ |l0[ |ll|, O. In particular, 
the Unified Approach proposed in Ref. |g] has been widely publicized by the Particle Data Group |13|] and used by 
several experimental collaborations. 

Using the Unified Approach we can derive confidence intervals for Ra- Figure [lj shows the confidence belts in the 
Unified Approach for a normal distribution with unit standard deviation for 90% (1.64a - ), 99% (2.58a), 99.73% (3a) 
and 99.89% (3.06a) CL. One can see that the measured value (|10| ) of R e A p /a A * p implies that 

0< -g^r < 6.32 at 99.89% CL, (11) 



i.e. active v^ or v T are present in the solar neutrino flux on Earth at 99.89% CL. Equation (O) implies that there 
is a 0.11% probability that the true value of Ra/ct^ 1S zer0 or l ar g er than 6.32. This probability is the same as the 
probability of a fluctuation larger than 3.06<r calculated in Ref. M assuming Ra = 0. However, our result have been 
derived without making any assumption on the true unknown value of Ra and has a well defined meaning in the 
framework of Frequentist Statistics: whatever the true value of Ra, the interval (O) belongs to a set of intervals that 
could be obtained in the same way from repeated measurements and have the property that 99.89% of these intervals 
cover the true value of RA/cr c A xp . 

In order to derive a lower limit for the averaged flux of v^ T on Earth, we consider in the following 99% confidence 
intervals. From Fig. [l] we obtain 

0.74 < 44 < 5.63 (99% CL) , (12) 

&A 

whose meaning is that there is a 99% probability that the interval ( |l2|) covers the true unknown value of Ra/ct a * p . 

For /b (Pi/ e _n/ r )j that gives the flux of active v^ T averaged over the common Super-Kamiokande and SNO response 
function, relative to the SSM 8 B v e flux, we find 

0.17 < / B (P„.-„„J < 1-26 (99% CL) . (13) 



Hence, we can say that the averaged flux of v^, T on Earth is larger than 0.17 times the 8 B v e flux predicted by the 
Standard Solar Model at 99% CL. This is an evidence in favor of relatively large v e — > v^ T transitions if fg is not 
too large. 

One could argue that it is possible to derive a more stringent lower limit for /b (P Vt ^ v t ) by calculating a confidence 
belt without left edge, instead of the one in Figure p] calculated in the Unified Approach. Such a procedure is not 
acceptable, because it would lead to undercoverage if not chosen a priori, independently from the data, as shown in 
Ref. p| for the case of upper limits. The correct procedure is to choose a priori a method like the Unified Approach 
that gives always sensible results and apply it to the data, as we have done here. A priori one could have chosen 
another method, as those presented in Refs. (lfl, O, O], that may have even better properties than the Unified 
Approach |lj, Q, but we have verified that the intervals (|rj])-(|l3|) do not change significantly. 

Unfortunately, we cannot derive a model independent lower limit for the averaged v e — ► i/„ T probability (P v _^„ T ), 
because /b could be large. However, from Figure [| we can say that Ra/o^ v > at 99.89% CL (see Eq. (pi])), and 
hence 



P„ e ^„_ >0 at 99.89% CL 



(14) 



in the range of neutrino energies covered by the common SNO and Super-Kamiokande response function presented in 
Ref. §. 

On the other hand, it is interesting to note that the relations (fil) and (0) allow to derive a model-independent lower 
limit for /b, taking into account that 



<P„.^„,., r )<l-<P„ a ^„.>- 
Using this inequality, from Eqs. (|l|) and Q2j) we obtain 

, ^ (gvJ p ( (°V e ) ,\n _, 

/B > -, r -flSK — -, r — i- -ftSNO = JB.min 

(<H,r) \K,x) / 

From Eqs. (||), (g) and (Q), the experimental value of /B,min is 

& = 1-031 ±0.197. 



(15) 



(16) 



(17) 



Since the central value of f^^m 1S 5-2c awav horn zero, we can calculate the resulting 99% CL interval for /B.min 
using the Central Intervals method (see jl3|), that gives the same result as the Unified Approach far from the physical 
boundary /B,min > 0. Since in the Central Intervals method 99% CL corresponds to 2.58<r, we obtain the confidence 
interval 

0.52 < JB, m in < 1-54 (99% CL) . (18) 

Therefore, we can conclude that the SNO and Super-Kamiokande data imply the model-independent lower limit 

/b>0.52 (99% CL). (19) 



This is a very interesting information for the physics of the Sun. 

So far we have not considered the possible existence of exotic mechanisms that produce v e 
transitions 
neutrinos 4 
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addition or alternative to v e — + v^ T transitions), such as resonant spin-flavor precession of Majorana 
In this case, Eq. (En must be replaced with 



-Rsk — /b {Pv, 



,/b 



(^.t) 



(Pv 



„r) + ^4 <?*.-«„) + 7?4 Pur**.) 



Wu e ) 



and Eq. (||) with 



R A = h 






(*».) 



(0-v.) 



(20) 



(21) 



3 We would like to thank a referee of the first version of this paper for pointing out the possibility of u e — ► v^ t r transitions of solar 
neutrinos. 



In the case of Majorana neutrinos the right-handed states are conventionally called antincutrinos. 



Using the 8 B neutrino spectrum given in Ref. p8[ , the neutrino-electron elastic scattering cross section calculated in 
Rcf. |H| taking into account radiative corrections, and the Super-Kamiokande energy resolution given in Ref. |2(J, we 
obtain the following values for the ratios of the averaged cross sections in the Super-Kamiokande experiment for the 
threshold energy T e SK = 8.60 MeV: 

y»'") =0.114, -^4 = 0.120. (22) 



Hence, we have the useful inequalities 



(ov 



i^,r) < pA < ^!±A, (23) 



The lower bound in Eq. dill) i mplies the existence of solar v e — > z^. r or z^ e — ► z^, jT or ^ e —* v e transitions at 99.89% 
CL. The inequalities in Eq. (|12[) imply that the quantity on the right-hand side of Eq. (|2l| ) is limited in the interval 
(0.025,0.19) at 99% CL. Using the inequalities Q, we obtain 



0.17 < / B [{Pu^u^) + (i^.-p MlT > + (-P.,.-*.}] < 1-67 (99% CL) . (24) 

Therefore, the averaged flux of v^, v T , v^, v T and v e on Earth is larger than 0.17 times the 8 B v e flux predicted by 
the BP2000 Standard Solar Model at 99% CL. 

Let us derive now the most general model-independent lower limit for /b (assuming only that the Super-Kamiokande 
and SNO events are produced by neutrinos or antineutrinos generated as v e from 8 B decay in the Sun). Using the 
inequality 

{Pu e ^, T ) + (Pv e ^, T ) + {Pu e ^u e ) < 1 - (Pu.-*u.) (25) 

and those in Eq. (|23|), from Eqs. ([]]) and ( |20| ) we obtain again the limit in Eq. (fl6|). Therefore, Eq. ( |19|) gives the 
most general model-independent lower limit for /b following from the SNO and Super-Kamiokande data. 

In conclusion, we have considered the model independent relations ([[]), (g) O, ||, [7| (and ([j]), (|20J) ) and the rates 
measured in the SNO m and Super-Kamiokande |2| solar neutrino experiment in the framework of Frequentist 
Statistics. We have shown that the probability of v e — > i^ )T (and z^ e — > 9^ T , v e — > j/ e ) transitions is larger than zero at 
99.89% CL in the range of neutrino energies covered by the common SNO and Super-Kamiokande response function. 
We have found that the flux of v^^ (and v^ T , 9 e ) on Earth averaged over the common SNO and Super-Kamiokande 
response functions is larger than 0.17 times the 8 B v e flux predicted by the BP2000 Standard Solar Model at 99% 
CL. We have derived a model-independent lower limit of 0.52 at 99% CL for the ratio /b of the 8 B v e flux produced 
in the Sun and its value in the BP2000 Standard Solar Model [§. 
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FIG. 1: Confidence belts in the Unified Approach |9[ for a normal distribution with unit standard deviation. The regions 
between the solid, long-dashed, dotted and dash-dotted lines correspond, respectively, to 90% (1.64a), 99% (2.58a), 99.73% 
(3a) and 99.89% (3.06a) CL. The thick solid vertical line represent the measured value of R^ p /a™ p (Eq. ([f0[)). 



